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FOREWORD 


This report presents results of the expansion and improvement of the 
FORMA system for response and load analysis. The acronym FORMA stands 
for FORTRAN Matrix Analysis. The study, performed from 16 May 1975 
through 17 May 1976 was conducted by the Analytical Mechanics Department, 
Martin Marietta Corporation, Denver Division, under the contract NAS8- 
31376. The program was administered by the National Aeronautics and 
Space Administration, George C. Marshall Space Flight Center, Huntsville, 
Alabama under the direction of Dr. John R. Admire, Structural Dynamics 
Division, Systems Dynamics Laboratory. 

This report is published in seven volumes: 

Volume I - Programming Manual, 

Volume IIA - Listings, Dense FORMA Subroutines, 

Volume IIB - Listings, Sparse FORMA Subroutines, 

Volume lie - Listings, Finite Element FORMA Subroutines, 

Volume IIIA - Explanations, Dense FORMA Subroutines, 

Volume IIIB - Explanations, Sparse FORMA Subroutines, and 
Volume me - Explanations, Finite Element FORMA Subroutines. 
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ABSTRACT 


This report presents techniques for the solution of structural 
dynamic systems on an electronic digital computer using FORMA (FORT RAN 
^totrix ^alysis). 

FORMA is a libiary of subroutines coded in FORTRAN IV for the effi- 
cient solution of structural dynamics problems. These subroutines are 
in the form of building blocks that can be put together to solve a large 
variety of structural dynamics problems. The obvious advantage of the 
building block approach is that programming and checkout time are limi- 
ted to that required for putting the blocks together in the proper order. 

The FORMA method has advantageous features such as: 

1. subroutines in the library have been used extensively for many 
years and as a result are well checked out and debugged; 

2. method will work on any computer with a FORTRAN IV compiler; 

3. incorporation of new subroutines is no problem; 

4. basic FORTRAN statements may be used to give extreme flexi- 
bility in writing a program. 


Two programming techniques are used in FORMA: dense and sparse. 
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LIST OF SYMBOLS 


matrix 

column matrix . 

\ vector 
row matrix ^ 

transpose (when symbol is a superscript) 

m designates the row size of matrix 
n designates the column size of matrix 

a designates an element of matrix [a] 

1 designates the ith row of matrix [aJ 
j designates the i th column of matrix [a] 
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!• INTRODUCTION 


This volume presents an explanation of the function of each sparse 
subroutine in the FORMA library. Example problems are given in some 
cases to clarify the operations performed by a subroutine. 
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II. SUBROUTINE EXPLANATIONS 

The subroutines are given in alphabetical order with numbers 
coming before letters. 



YAA 


Subroutine YAA calculates the multiplication of a matrix, In FORMA 
sparse notation, by a scalar. In matrix notation, 

NRxNC “ “ NRxNC 


where 



NR is the number of rows of each matrix, and 
NC is the number of columns of each matrix. 


EXAMPLE 


Consider input of « “ 2; and 



The reader can easily verify the output to be 



DESCRIPTION OF TECHNIQUE 

[a 3 and [z ] are stored on utility tapes NUTA and NUTZ as row partitioned 
matrices in FORMA sparse notation. 

Matrix [A] is read from NUTA and each non-zero element is multiplied by a. 
Each partition of [A 3 then becomes a partition of [z] and is stored on 

Nirrz. 



YAABB 


Subroutine YAABB calculates the summation of two matrices, each 
matrix multiplied by a scalar. Both matrices are In FORMA sparse 
notation. In matrix notation, 




where 


2 “ <x a 

ij iJ 


+ 6 b 


iJ 


1 = 1, NR 
j = 1. NC 


NR Is the number of rows of each matrix, and 
NC Is the number of columns of each matrix. 


The number of rows of [A] and [B] must be equal, and the number 
of columns of [A] and [B] must be equal. 

DESCRIPTION OF TECHNIQUE 

lA], [B], and [Z] are stored on utility tapes NUTA, NUTB, and NUTZ 
as row partitioned matrices In sparse notation. 

Partitions of the matrices are operated on In core In the V and 
LV work spaces as follows: (1, m, n refer to partition number) 


LV 


[A] a) 

[B] (m) 
[2] (n) 


1=1, NPARTA 
m = 1, NPARTB 
n = 1, NPARTZ 


The locations of elements In [A] and [B] are compared with each 
other. If the locations are equal, elements a and b are summed 
and stored in (Z]. If the locations are not equal, the lesser 
location element is stored in [Z]. 



YASSEM - 1/2 


Subroutine YASSEM places (assembles) a matrix [A] into a se< ond matrix [ 2 ] 
starting at a designated row, column location (IRZ, JC2, respectively) in 
[zj. Both matrices are in FORMA sparse notation. The elements of [aJ will 
replace corresponding elements in the original [z] . Before the first use 
of this subroutine in forming [z] , it is important that [z] is correctly 
defined. For example, if [zj is to be originally null, subroutine YZERO 
should be used. This subroutine may be called repeatedly to form [z] from 
the assembly of several Ta] matrices. The [A] matrix must be within the 
row, column limits of fzj . In subscript notation, 



where 

i » k + IRZ -1; 
j - + JCZ -1; 

NRA is the number of rows in [Al ; and 
NCA is the number of columns in [A] . 


EXAMPLE 


*3x4 


[z] 

Matrix [a] 


matrix 

defined as 




1. 0. 

0. 0. 


0. 2. 

0. 0. 

* 

0. 0. 

0. 0. 


- 

- 



3. 4. 5. 

2x3 “ 

6. 7. 8. 


is to be assembled 


into [z] starting at the 2, 1 location (IRZ * 2, JCZ = 1) of [Z] 
The result of this operation will be 


[z] 


3x4 


a 


1 . 0 . 0 . 0 . 
3. 4. 5. 0. 

6. 7. 8. 0. 


DEXRJ.PTION OF TECHNIQUE 

[aJ and [Z^ are stored on utility tapes NUTA and NUTZ as row partitioned 
matrices in FORMA sparse notation. 
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Matrix [a] is read from NUTA one partit;..on at a time. The location 
numbers of its non-zero elements are revised by adding the quantity 
[100000* (I RZ- 1 ) + (JCZ-l)l t;o them and the revised matr^A [aJ is stored 
on NUTl. 

In the next operation, matrix [z 3 is read from NUTZ one partition at a 
tln^. Any non-zero elements existing in the area of the matrix where 
[AJ is to be inserted, are eliminated and [zj is stored on NUTZ. 

Matri;. KJ is read from NUT 2 and written on NUTZ. Matrix [Xj is then 
read from NUTl and stored with [zj on NUTZ. 

The final step is a_complished by calling subroutine YLORD to correctly 
order the remaining non- zero elements of (zJ together with the revised 
elements of [A] , and the resulting matrix is stored on NUTZ. 



YBSL3A - 1/2 

Subroutine YBSL3A solves the linear simultaneous algebraic equa- 
tions [A] [Z] a [B], where [A] has been previously decomposed 

Into tU)^ fDj tU]. In matrix notation. 


Decomposition of [A] gives 

[Uj*^ fDj [U] [Z] « [B] [2] 

Letting 

[U] [Z] = [Y] [3] 

and 

tDj [Y] = [G] (4] 

gives 

lu]'^ [G] = [B]. [5] 

Because 


Because 



the eleuents of [G] can be solved as 



k=l 
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Using [4] we can now solve for [Y] as 





Now, because 




^ "12 

^3 

• • • “in 


1 

“23 

• • • 

[U] - 


1 

• * • “3N 


0 

1 


Che back solution for [Z] Is 



k=l+l 


Subroutine YBSL3A uses [U] and [OJ Input, but expects them In a 
special format. They are not In FORMA sparse notation. Input 
matrix [B] and output matrix [Z] are In standard FORMA sparse 
notation. 



YBTAB - 1/2 


Subroutine YBTAB calculates the triple matrix product for matrices in 
FORMA sparse notation. In matrix notation 

“ f^^NCBxNRB ^^^NRBxNRB f®^N8BxNCB 

where 

"•ij “ Hi Hk \j (j = l\ NCb) • 

MRB is the number of rows of [B^ and the size of [a 3 (square) . 

NCB is the number of columns of [bJ and the size of [zj (square) . 

Because the number of columns of [A^ must be eaual to the number of rows 
of [bT , and because the number of columns of f^BjJ (or the number of rows 
of L B j ) must be equal to the number of rows of ^3 > then [a 3 must be a 
square matrix. The answer [z3 is likewise a square matrix. 

Theorm ; If [a 3 is syranetric (that is, [a 3 * W^). tl»en [z3 is symme- 
tric. 

Proof : [z3^ = ([sf [a3 [b] 

. [b]I [a]^ ([bB^ 

= [b]^ [a] [b] 

' [B] . 

DESCRIPTION OF TECHNIQUE 

[a3 . [b 3 , and [z] are stored on utility tapes NUTA, NUTB, and NUT2 
respectively as row partitioned matrices in FORMA sparse notation. 

Partitions of matrix [b] are read from NUTB and the location numbers 
of the non-zero elements are transposed. These partitions now are 
partitions of Tb]! and are stored on NUTZ. 

Subroutine YMULT is called to form the product of [b]^ [a] . This 
product is stored on NUTl. 
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If (a] l8 not symmetriCt subroutine YMULT is called to complete the 
triple matrix product |a]) (B] , and the resulting [Z3 is stored 

on NUTZ. 


If Is symmetric, the result of [B] [A] is read into core, one par- 
tition at a time, from NUTl. Each row of CAJIs postmultiplled by 

giving a row of QzJ . To save computer time, only the for j ® l,i 
are formed giving the lower half of the symmetric matrix [ZJ , QzJ is 
stored on NUTZ. The logic used in this portion of YBTAB is a slighu modi- 


fication of the logic used in YMULT. 



YDCM3A - 1/4 

Subrou*’ins 1.DCM3A decomposes [A] a matrix in FORMA sparse nota- 
tion to f > m an upper triangular matrix with ones on the diagonal 

[U] and a diagonal matrix fDj such that [A] * (Uj^ fDj [U]. [A] 

must br reil, square, and symmetric. The decomposition technique 
is attributed to Gauss. In matrix notation 



N is the size Tf the matrices (square): 

DESCRIPTIC ^ - OF TECHNIQUE 

To de^ ..;lne the elements of [U] and fDj , first consider the 
mult ‘plication of two latrices. 


[A] = [c] ru] 

where [C] is a lower triangular matrix of the form 


( 2 ) 



YDCM3A - 2/4 



and [U] has been previously defined. 
Because [A] Is symmetric, 

[A] » (A]"^ 


=■ [uf tc]*^ 



Due CO Che uniqueness of the factorization of [A] Into the product 
of the two triangular matrices, one of which has ones on the diagonal. 
It follows chat 
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By defining 


m - 


‘=11 

0 

... 0 

0 

‘=22 

... 0 

0 

• O 

... c, 

- 




NN 


Equation [1] Is thus derived from Equation [2]. 

By performing the matrix products Indicated by Equation [1], the 
elements in [A] are obtained generally as 

“ij ■ “ll '*11 “ij * “21 ‘*22 “2J ‘‘ii “ij 

and 

“ii " ‘‘a “u '^22 “a + ••• ■" '*11 

From these two equations we calculate 
‘^ii ” ^11 




1-1 

“'ll ' “ii ■ S '^1 

k=l 
1-1 


(J = 2. N) 
(1 = 2, N) 


u,. - |a_ - ^ u,^, d,,,, (i < j) 


'ij \"ij ki -kk “kj/fii 

k-1 


Uii o 1.0 


Ufj = 


(1 = 1, N) 
(i > J) 


Because of the shape of [A] and the corresponding shape and 
density of [U], this decomposition is performed using banded 
programming logic. [U] and not stored in FORMA sparse 

notation but are in a form to be used as input to subroutine 
YBSL3A. 
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MISCELLANEOUS 


The diagonal elements of TDJ are the determinant ratios of [A]. 

The determirant of [A] Is the product of these determinant ratios. 
That la, 

N 

|a1 “ TT 

i=l 

REFERENCE 


V. N. Faddeeve, Computational Methods of Linear Algebra. Dover 
Publications, Inc., New York, 1959. 



YDC0M2 - 1/4 

Subroutine YDC0M2 decomposes (A] a matrix In FOKMA sparse nota- 
tion to form an upper triangular matrix with ones on the diagonal 

[U] and a diagonal matrix fDj sucli that [A] = (U]^ [\)] (U). (A) 

must be real* square, and symmetric. The decomposition technique 
Is attributed to Gauss. In matrix notation 




where 



N is the size of the matrices (square): 

DESCRIPTION OF TECHNIQUE 

To determine the elements of [U] and fDj , first consider the 
multiplication of two matrices. 


lA] = [Cl [U] 

where (C) is a lower triangular matrix of the form 


( 2 ) 



YDCQM2 - 2/4 



aud [U] has been previously defined. 
Because [A] is symmetric, 

[A] = lA]"^ 

= lU]^' [CJ^ 



Due to the uni cjueness of the factorize t ion of [Aj into tlie product 
of the two triangular matrices, one of which has ones on the diagonal, 
it follows that 
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By defining 


[DJ = 


‘^ll 

0 

o 

0 

"22 ••• 

0 

0 

0 

"nn 


Equation (IJ Is thus derived from Equation [2]. 

By performing the matrix products indicated by Hcjuation [IJ, 
elements in (A] are obtained generally as 

a_ u, . d- , u. . + u^ . d^^ . + ... + cl . . u . . 

ij ii 11 1) 2i 22 2j 11 ij 

and 

2 2 

a = d. - u- . + d^.. . + ... + d., 

li 11 li 22 2i li 

From these two equations we calculate 
"*11 " “'ll 


“ij ■ 


1-1 


dii = a 


u-E 


, 2 

! Al- 
lele ki 


k-1 

i-l 


u , , “ 1 a . . - 7 u 

ij \ 1] Z-/ 


ki “kk \j rn 


(j = 2. N) 

(i = 2, N) 

(i ^ .1) 


k-1 


li ^ 

(i = 1. N) 

, . =0.0 

(1 > j) 

Ij 


Because of the shape of [A] and the corresponding shape and 
density of (U], this decomposition is performed using banded 
programming logic. 


the 
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For the banded programming logic employed In YDC0M2, only half of QAJ is 
stored. Rows of are stored \dilch Include only the diagonal element 
through the last non-zero in the row. Output matrices Qu] and are 
written on NUTU and NUTD In FORMA sparse notation. 


MISCELLANEOUS 


The diagonal elements of determinant ratios of QA3 . The 

determinant of C^3 product of these determinant ratios. That Is, 



N 

TT 


1=1 


d 


11 


REFERENCE 

V. N. Faddeeve, "Computational Methods of Linear Algebra". Dover 
Publications, Inc., New York, 1959. 
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Subroutine YDC0M3 decompose? [A] a matrix In FOKMA sparse nota- 
tion to form an upper triangular matrix with ■'nes on the diagonal 

[U] and a diagonal matrix such that [A] = (UJ^ fl)J [U]. |A) 

must be real, square, and symmetric. The decomposition technique 
if? attributed to Gauss. In matrix notation 



N is the size of the matrices (square): 

DES CRIPTION OK TECHNIQ UE 

To determine the elements of [U] and fDj , first consider the 
multiplication of two matrices. 


lA] = [C] [U] 

where (C) is a lower triangular matrix of the form 


12 ) 




and [U] has been previously defined, 
because [A] is symmetric, 

[AJ = [A]"^ 

= [u]^' [:]^ 



Due. to the uniqueness of the factorization of [Aj into the product 
of the two triangular matrices, one of which has ones on the diagonal, 
it follows that 
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By defining 


fDj 


c 0 
11 


... 0 


^-22 ••• ^ 


'NN 


Equation (1] is thus derived from Equation (2). 

By performing tlie matrix products indicated by Equation lll» the 
elements in fAJ are obtained generally as 

a. . * u, . d, , u, . + u«, d-- u., . + . . . + d. , u. . 
ij li 11 Ij 2i 22 2j 11 ij 

and 

®ii “ *^11 *^22 “a ^ii 


From these two equations we calculate 


d 


11 


a 


11 


u 


IJ 


a ./d„ 

Ij 11 


(J = 2, N) 



Because of the shape of [A] and the corresponding shape and 
density of fUj, this decomposition is performed using banded 
programming logic. 
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For the banded programnlng logic eaployed in YDC0M3, only half of CA3 i> 
stored. Colunna of [a 3 are stored which include only the top non- zero 
down to the diagonal element of the column. Output matrices C^l and 
are written on NlTfU and NUTD in FORMA sparse notation. 


MISCELLANEOUS 


The diagonal elements of are the determinant ratios of ([Aj . The 
determinant of C^] is the product of these determinant ratios. That is. 



REFERENCE 

V. M. Faddeeve, Conq>utational Methods of Liiwar Algebra. Dover Publi- 
cations, Inc., Mew York, 1959. 
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Subroutine YMAG places the eleiaents from a vector (row or column matrix 

1 A I into the corresponding diagonal locations of a square matrix • 

A I and are stored on NUIA and NUIZ in FORMA sparse notation. In 
subscript notation, 


z 


ii 




z 


ij 


» 0 . 


(i » 1, N) 
(i ^ j) 


The ate not stored. 
In matrix notation 


[ 2 ] 


NxN 




where H is the size of (square), and the length of {a| . 


EXAMPLE 

Consider input of 
where N ■ 3. 

The result of this subroutine will give 

1 . 0 . 0 . 

0 • 2 • 0 . • 

0. 0. 3. 




DESCRlPnON OF TECHNIQHE 

Vector |aJ is read from NUTA, one partition at a time. The location 
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nuabers for the non*zero elements of |a| are modified to reflect diagonal 
elements. These partitions are then mitten on NUTZ to form . 



YDISA - 1/2 


Subroutine YDISA removes (disassembles) a matrix from matrix ^A) 
starting at a designated rov, column location (IRA, JCA, respective!) 
In [a] . The Qz] matrix must be within the row, column limits of CaO . 
Matrices [A] and Qz] are stored in FORMA sparse nat tlon. In subscript 
notation, 

/i = I, NRZ\ 

\j = 1. NCZ j 

k « 1 + IRA - 1 
i = j + JCA - 1 

NRZ is the number of rows of Cz} , and 
NCZ is the number of columns of Cz3 . 

EXAMPLE 

Consider a matrix defined as 

<■ 

1 . 0 . 0 . 0 . 

3. 4. 5. 0, 

6 . 7 . 8 . 0 . 


Matrix ^^^2x3 obtained from CaJ starting at the 2,1 location 

(IRA = 2, JCA = 1) of CaJ • The result of this operation will be 



The matrix [A] remains as originally defined. 
DESCRIPTION OF TECHNIQUE 



"ij “ \i 

where 


Partitions of [a] are read from NUTA into the first quarters of the V 
and LV workspaces. The element locations of [A] are searched for 
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locations In the range of • Any elements found in this range are 
moved to the second quarter of V and location numbers for these elements 
are formed In the second quarter of LV. The data in the second quarters 
of V and LV become partitions of vhlch Is stored on NUXZ. 



YDTOS 


Subroutine YSTOD converts a dense matrix In core to a FORMA sparse 
matrix written on a peripheral device. Only the nonzero elements 


of dense matrix {A] are stored in sparse matrix [Z]. 


For each a 


IJ’ 


a location number is formed as 


LV(k) » 100000 * i + j and 


V(k) «= a^j. 

LV(k) and V(k) for k <= 1, NNZA are then written on the peripheral 
device, where NNZA Is equal to the number of nonzeroes In [A]. 



YIN 


Subroutine YIN is used to transfer real data from a computer 
peripheral device Into core. Utilization of this subroutine enables 
a programmer to easily take advantage of computer -dependent input 
routines which may be more efficient than FORTRAN read statements. 
Data transmitted using subroutine YOUT may be retrieved using sub- 
routine YIN. 



YINI 


Subroutine YINI is used to transfer integer data from a computer 
peripheral device into core. Utilization of this subroutine en- 
ables a programmer to easily take advantage of computer dependent 
input routines which may be more efficient than FORTRAN read state- 
ments. Data transmitted using subroutine YOUTl may be retrieved 
using subroutine YINI. 



YLORD 


The logic in FORMA sparse subroutines depends on the fact that 
the locations of elements In the matrices are In order by row and 
column. 

Subroutine YLORD orders all nonzero matrix element locations and 
the elements. Ordering Is defined as follows: all elements of 

a row are ordered by column from left to right; all rows are then 
ordered from the first through the last. 

DESCRIPTION OF TECHNIQUE 

Because no partition of a FORMA sparse matrix exceeds KV/4 In 
length, four partitions of [A] are brought into core together and 
the element locations and elements of all four are ordered simul- 
taneously. The method of ordering employed in this process is 
derived from the method of R. C. Singleton.* All partitions of 
[A] are ordered, four at a time, using this process. The ordered 
partitions then are merged together to form the completely ordered 
matrix [Z]. 


*R. C. Singleton: An Efficient Algorithm for Sorting with Minimal 

STorage. Research Memorandum, SRI Project 387531-132, September 
1968. 



■5fMQD2A 


Subroutine 'XM0B2A calculates the mode shapes and natural fre- 
quencies of a structural model described as 

K] [♦] - [M] [^] 11] 

where 

[KJ » [Kj - Ag [Ml 
and 

. Ag [I]. 

If Ag - 0, then [K] = [K] and and [1] becomes the 

more familiar eigen problem 

tK] U] * [M] [$] Tu)2J. 

Subroutine 'SMOD2A uses the composite structure iterative Rayleigh- 
Rltz method of Dr. John Admire. 

If initial Raylelgh-Ritz coordinate displacements [qj are avail- 
able, they may be input to the subroutine. Otherwise random numbers 
will be used. 

Input matrices [K], [M], end [qj and output matrix [(frj are in 
FORMA sparse notation. 
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Subroutine !2MULI calculates the product of two matrices stored in FORMA 


sparse 

notation. 

'‘.n matrix notation, 

^•^^NRAxNCB " 

^^^NRAxNRB 

f®^NRBxNCB 

where 






NRB 

E 

k«l 

®ik 


/i - 1, NBA 
1 j = 1, NCB 


NBA is the number of rows of fA] and [Z^ , NRB is the number of rows of 
QbJ and the number of columns of Ca 3 • NCB is the number of columns of 
£B] and [zj . The number of columns of must be equal to the number 
of rows of Cb3 • 

Theorem ; Multiplication of matrices is not commutative in general. 

That is. 


[A] W t [B] [A] 


for any values of and b^^j 

Theorem ; Multiplication of matrices is associative. That is, 

[a] ([b] [c]) = ([a][b]) [c]. 

Theorem ; Multiplication of matrices is distributive. That is, 
[A] ([B] + [C]) » [a] [b] + [A] [C 3 . 

EXAMPLE 

Consider input of * 

The reader can easily verify the output to be 



- [27. 14. 33.] . 
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7JESCR1PTI0N OF TECHNIQUE 

Partitions of are read from NUTA Into the first quarters of the V 
and LV workspaces. Partitions of [Bj are read into the second quarters 
of the workspaces and are premultiplied by i row of [A] . A dense row 
of Cz] is formed in the third quarter of V based on the row of [Aj mul- 
tiplied by all of matri-x [BJ . When the dense row of [Z] is completely 
formed, the non-zeros are transferred to the fourth quarter of V and an 
appropriate vector of location numbers is formed in the fourth quarter 
of LV. The data generated in the fourth quarters of V and LV become 
partitions of [zj and are written on NUTZ. 

Limitation: 

One fourth of the dimension size of V and LV, (K\/4), must exceed or be 
equal to NCB and NRB. 



YMULTl - 1/2 

Subroutine YMULTl calculates the special product of two FORMA 
sparse notation matrices. In matrix notation, 

^^^NRAxNCB " ^^^NRAxNRA ^®^NRAxNCB 
NRA 

‘ij " S “ik ‘■kJ 

k=l 

NRA is the number of rows of [A], [B], and (Z] where I'A) is 
square. NCB is the number of columns of [B] and [Z]. NRA is 
assumed to be greater than NCB. 

[B] and therefore [Z] are assumed to have no elements equal to zero. 

[A], [B], and [Z] are scored on utility tapes as row partitioned 
matrices In sparse notation. 

Partitions of the matrices are operated on in central memory as 
follows (I, m, and n refer to partition number): 

V LV 


/l = 1, nra\ 
\j = 1, ncb/ 


KV/4 — — 


I^A(i) i = 1, NPARTA 

KV/2 — . 


B(m) m =« 1, NPARTB 

not used 
KV 


Z(n) n = 1, NPARTZ 


not used 


where V and LV are dimensioned work spaces and KV is the dimension 
size of V and LV. 


The number of tvonzero terms in each partition of (Z), NNZPZ, is 

> I 


f tnj\ 

^1/ncb * ncb. 


The number of rows in each partition of [Z], NRPZ, is NNZPZ/NCB. 

The number of partitions of [Z], NPARTZ (-NPARTB) , is (NRA-1)/ 
NRPZ + 1. 
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^ — 




— 

A, 




2 , 

1 




1 



B 

= 

» * * • 

A 




• * • « 

z 

nparta 

— 




npartz 


NFRPZ(l) » 1 
NLRPZ(l) 

NFRPZ(^IPARTZ) 
NLRPZ(NPARTZ) = NRA 








i = NFRPZ(fc),NLRPZ(?.) 
J = 1, NCB 
k = 1, NRA 
9. = 1, NPARTZ 


Where 


NFRPZ is the number of the first row in a partition of Z, and 
LLRPZ is the number of the last row in a parti cion of Z. 



YMULT2 


Subroutine YMULT2 calcilates the special product of two FORMA 
sparse notation matrices, in matrix notation, 

T 

^^^NCBxNCB ^ ^^^NCBxNRB ^^^NRB^NCB 
NRB 

k-1 

NRB is the number of rows of [A] and [BJ and the number of columns 

T T 

of [A) • NCB is the number of rows of [A] and [Z] and the number 

of columns of [A], [B], and [Z]. [A], [Bj, and [Z] are assumed to 

have no elements equal to zero and [Z] is symmetric. 

Statements from FORMA subroutine ATXBB2 are used in this subroutine* 


/i = 1, NCB\ 
\j = 1, NCB/ 



YMULT4 


Subroutine YMULTA calculates the special product oT a FORMA sparse 
notation matrix postmultlplled by a square matrix stored in core. 
In matrix notation » 

^^^NRAxNCB “ ^^^NRAxNCB ^®^NCBxMCB 


z 


ij 


NC3 


E. 


*ik 


/i = 1, NRA\ 

\j = 1, ncb/ , 


where 

NRA Is the nunber of rows of (A] and [Z], and NCB is the number of rows 
of [B] and the number of columns of [A], [B] and (Z]. NRA is greater 
than NCB. 

(A] and [Z] are assumed to have no elonents equal to zero. They 
are stored on utility tapes in FORMA sparse notation. 
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Subroutine IMULTA calculates the product of two matrices stored in FORMA 
sparse notation* In matrix notation 


f^^NRAxNCB 

II 

^®3nrBxMCB 

where 




z. . = 

MSB 

E 

k*l 

1 • 
ik kj 

/i = 1, NBA 
\i = 1, NCB 


NRA is the number of rows of Ca] and [z] . NRB is the number of rows 
of the number of columns of QaJ « NCB is the number of columns 

of ([BJ and [Zj . The number of columns of ^a 3 must be equal to the 
number of rows of [BJ « 

Theorem ; Multiplication of matrices is not comnutative in general. 

That is, 

[a] [b] f [b] [a] 


Theorem: 


for any values of a., and b... 

^ ij ij 

Multiplication of matrices is associative. That is. 


[A] (W tc]) = ([A] W) [C] . 


Theorem : 


Multiplication of matrices is distributive. That is. 


[a] ([b] + [c]) = [a] [b] + [a] [c] . 


EXAMPLE 

Consider input of [a]ij^ 2 * 2*3 ^2x3 ' 

The reader can easily verify the output to be 


[Z] 


1x3 



- 8 . 

11 . 


33.] 




- 8 . 

11 . 


9. 

12 . 


] 


[ 1 . 

[27 


14. 
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DESCRimON OF TECHNIQUE 

Partitions of [aJ are read from NUTA into the first quarters of the V 
and LV workspaces. Partitions of are read into the second quarter 
of the workspaces and are pretnul tip lied by a row of [A] . A dense row 
of formed in the third quarter of V based on the row of fA] mul- 

tlplied by all of matrix [bJ . When the dense row of fZj is completely 
formed, the non- zeros are transferred to the fourth quarter of V and an 
appropriate vector of location numbers is formed in the fourth quarter 
of LV. The data generated in the fourth quarters of V and LV become 
partitions of [Zj and are written on NUTAZ. Matrix fA3 is, therefore, 
replaced by matrix (^Zj in this subroutine. 

Limi tation: 

One fourth of the dimension size of V and LV, (KV/4), must exceed or be 
equal to MCA and NCB. 
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Subroutine YMULTB calculates the product of two matrices stored in FORHA 
sparse not^.ion« In matrix notation 


^^^NRAxHCB “ 

t'^^NRAxNRB 

^^®^NRBxNCB 

where 





nbb 

z 

k=l 

3^, b, . 
ik kj 

II It 


NRA is the number of rows of and [Zj . NRB is the number of rows of 
fBj and the number of columns of . NCB is the number of columns of 

[B] and fzj . The number of columns of [A] must be equal to the number 
of rows of C®3 • 

Theorem ; Multiplication of matrices is not commutative in general. 

That is, 

[a] [b] # [b] [a] 

for any values of a,, and b,,. 

ij ij 

Theorem : Multiplication of matrices is associative. That is, 

[a] ([b] [c]) = ([a] [b]) [c] . 

Theorem : Multiplication of matrices is distributive. That is, 

[a] ([b] + [c]) - [a] [b] + [a] [c] . 

EXAMPLE 

Consider input of [A] 1x2 “ 2.] and ^2x3 = 

The reader can easily verify the output to be 



» [27. 14. 33.] 
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DESCRIPTIOK OF TECHNIQUE 


Partitions of CA] are read from NUTA Into the first quarters of the V 
and LV workspaces* Partitions of [Bj are read Ifito the second quarters 
of the workspaces and are premultiplied by a row of [A] • A dense row 
of Cz] Is formed in the third quarter of V based on the row of fAj mul- 
tiplied by all of matrix [,BJ . When the dense row of C^] is completely 
formed^ the non-zeros are transferred to the fourth quarter of V and an 
appropriate vector of location numbers is formed in the fourth quarter 
of LV* ^e data generated in the fourth quarters of V and LV become 
partitions of [Zj and are written on NUTBZ. Matrix fBj is therefore 
replaced by matrix Cz] in this subroutine* 

Limitation: 

One fourth of the dimension size of V and LV, (KV/4) , must exceed or be 
equal to NCA and NCB. 



YNOZEk 


Subroutine YNOZER scans each partition of a FORMA sparse matrix 
and deletes any terms equal to zero. The logic of sparse matrix 
subroutines becomes Inefficient If matrix terms equal to zero are 
stored and operated on. If prior matrix operations have computed 


values In [A] such that any a 
deleted. 


ij 


are equal to zero, 


these a 


iJ 


arc 



YOUT 


Subroutine YOUT la used to transfer real data from computer central 
memory to a peripheral device. Utilization of this subroutine 
enables a programmer to easily take advantage of computer-dependent 
output routines which may be mere efficient than FORTRAN write 
statements. Data transmitted using subroutine YOUT may be re- 
trieved using subroutine YIN. 



YOUTI 


Subroutine YOUTI is used to transfer integer data from computer 
central memory to a peripheral device. Utilization of this sub- 
routine enables a programmer to easily take advantage of computer 
dependent output routines which may be more efficient than FORTRAN 
write statements. Data transmitted using subroutine YOUTI may be 
retrieved using subroutine YINI. 



YPART 


Subroutine YPART partitions newly generated matrices and reparti- 
tions existing matrices in which the partition length is not 
compatible with the work -pace dimension (KV) of the current 
computer program. Partitioning ensures that no matrix row spans 
two partitions and the length of all partitions is less than or 
equal to one fourth of KV. 



YPNCHO 


Subroutine YPNCHO punches a matrix [Al onto cards using octal format to 
eliminate round-off error. Matrix [A] is stored on NUTA in FORMA sparse 
notation. Octal representation of the matrix elements is used because 
it gives an exact replica of the binary number used by a digital computer. 
A decimal representation will not give an exact replica. The matrix on 
punched cards is to be used only as an emergency backup for the matrix 
written on a storage tape. 

This matrix on cards is compatible with the input form for Subroutine 
YRSADO* A group of up to three consecutive elements from a row of the 
matrix are punched on each card. If all of the elements of a group are 
zero 9 punching of this card is suppressed. 

The first card punched contains the matrix name (in card columns 1-6), 
the matrix row size (in card columns 7-10), the matrix column size (in 
card columns 11-15), and the matrix shape (in card columns 16-21). This 
is followed by the matrix data. On any card of the matrix data, the first 
integer number (card columns 1-5) is the row number of the matrix elements 
on that card. The second Integer number (card columns 6-10) is the column 
number of the matrix element in the first data field (card columns 14-25). 
The next group contains octal numbers (up to three numbers in card columns 
14-25, 29-40, 44-55) that are the values of the matrix elements. This 
group of matrix elements is given in consecutive column order. The last 
card punched contains ten zeroes in card columns 1-10 to indicate the end 
of the matrix. 



YPUNCH 


Subroutine YPUNCH punches a matrix [a 3 of real numbers (a Fortran term 
for numbers with a decimal point) onto cards. Matrix [aJ Is stored on 
NUTA In FORMA sparse notation. This matrix on cards Is compatible with 
the Input form for Subroutine YREAO. A group of up to four consecutive 
elements from a row of the matrix are punched on each card. If all of 
the elements of a group are zero, punching of this card Is suppressed. 

The first card punched contains the matrix name (In card columns 1-6), 
the matrix row size (In card columns 7-10), the matrix column size (In 
card coluims 11-15), and the matrix shape (In card columns 16-21). This 
Is followed by the matrix data. On any card of matrix data, the first 
Integer number (card columns 1-5) Is the row number of the matrix elements 
on that card. The second Integer number (card columns 6-10) Is the colunn 
number of the matrix element in the first data field (card columns 11-27). 
The next group contains real numbers (up to four numbers in card columns 
11-27, 28-44, 45-61, 62-78) that are the values of the matrix elements. 
This group of matrix elements is given in consecutive column order. The 
last card punched contains ten zeroes in card columns 1-10 to Indicate 
the end of the matrix. 
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Subroutine YREAD reads a matrix of real numbers (a FORTRAN term 
for numbers with a decimal point) from either cards or tape and 
stores the matrix on a utility tape in FORi-a sparse notation. The 
matrix is then printed so that these input data are recorded with 
the arnwers of a run. A print suppression option iw available 
for a matrix read from tape. On option, the matrix read from 
either cards or tape may be written on a tape (by Subroutine 
YWTAPE) . 

The first data card read by Subroutine YREAD contains informatior 
to indicate whether cards or tape will be used. The information 
entered on this card (and subsequent curds for card input) is given 
below. The format is i identical to that used by Subroutine READ. 

Card Data Input Form 


Required entries are denoted by an * symbol below. Any other 
entry is optional 


Card Format 

Columns Type (1) Entry 


First Card 1-6 

7-10 

11-15 

16-69 

16-20 



Ou 

0 

0) 

D. 

H 

1 

0 ) 

4J 


72 


73-78 


t-i 


73-78 


A 

I 

I 

A 

A 


or 


A 


or 


^Matrix Name. Will appear in 
printout . 

^Matrix Row Size. 

^Matrix Column Size. 

Any remarks to further identify 
the input matrix. 

The word upper (or lower) must 
be specified if only the upper 
(or lower) half of a symmetric 
matrix is on the punched cards. 

$- Only if the Write-Tape is 
to be initialized by Subroutine 
INTAPE. The Write-Tape identi- 
fication will be from card 
columns 73-78. 

Anything other than $ is the 
Write-Tape is not to be ini- 
tialized. 

The Write-Tape Identification, 
(eg., T1234). Use with $ in 
card coluirn 72. 

REWIND. The Write-Tape will 
be rewound before being used. 
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(coAfdl) 


Card 

Coluims 

73-76 


73-78 

79-80 


Middle Cards 1-5 

6-10 
11-27 
28-44 
45-61 
62-78 

Last Card 1 -10 


Format 

Type (1) Entry 


or 


or 

I 

or 


LIST. The Write-Tape will be 
listed by Subroutine LTAPE 
after the matrix has been 
written on the Write-Tape. 
Anything else will be ignored. 
The Write-Tape Number (e.g., 21). 
Blank if the matrix is not to 
be written on tape. 


I 

I 

E 

E 


E 


E 

I 


*Row Nunijer of matrix elements 
on card. 

*Column Number of matrix ele- 
ment in first data field. 

*First data field with matrix 
elements. (2) 

^Second data field with matrix 
elements. (2) 

*Third data field with matrix 
elements. (2) 

*Fourth data field with matrix 
elements. (2) 

*Ten zeroes. 


Note (1) Format Type A allows any keypunch symbol. 

Format Type I allows only integer numbers right justified 
in the field. Format Type E allows only real numbers 
(a FORTRAN term for numbers with a decimal point) any- 
where in the field. 


Note (2) Only nonzero elements need be entered. 


Ari an example of 

card 

input 

to Subroutine YREAD consider the fol 

lowing matrix 

1 









0. 

3. 

0. 

6. 

sT 



0. 

2. 

4. 

0. 

0. 

0. 



0. 

7 . 

0. 

0. 

0. 

0. 

. 


This matrix is also to be written on tape number 21 that is to 
be initialized and identified as T4334. Figure 1 demonstrates 
how this information could be written on a coding form to facili- 
tate keypunching to cards. 



FORMA I 

pAe,f or 
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Figure 1 Example of Card Input for Subroutine YREAD 
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Tape Data Input Form 

Required entries are denoted with an * symbol below. Any other 
entry is optional. Only one card is used for each matrix read. 

Card Format 

Columns Type (1) Entry 


One Card 1-7 A 

10 


7-10 I or 


11-15 I 


16-21 A 


22-25 or 

21-27 or 

28-69 


^Narne of matrix to be read from 
the Read-Tape. 

Zero. The Read-Tape will move 
forward from its present posi- 
tion and search to the end of 
the tape. If the matrix is 
not found upon the first end- 
of-tape encounter, the tape 
will automatically rewind and 
make one more pass. If it is 
not found on the second end-of- 
tape encounter, an error message 
will be printed and the program 
will stop. 

Minus the location number of 
matrix on the Read-Tape. Tape 
vjill be positioned at the be- 
ginning of the location specified 
and then continue as described 
above for a zero in column 10. 

*The Read-Tape Number (e.g., 11). 
If positive, the matrix read will 
be printed in the output. If 
negative, the matrix read v;ill 
not be printed in the output - 

*Run number ot matrix to be read 
from the Read^Tape. 

REWIND, The Read-Tape will b - 
rewound bef‘'.e being used. 

LIST. The Read-Tape will be 
listed by Subroutine LTAPE. 
Anything else will be considered 
as part of the remarks described 
below . 

Any remarks to further identify 
the input matrix. 
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Format 

Type (1) Entry 

$* Only if the Write-Tape is 
to be initialized by Sub- 
routine INTAPE. The Write- 
Tape identification will be 
from card columns 73-78. 
or Anything other than $ if the 

Write-Tape is not to be ini- 
tialized. 

A The Write-Tape identification 

(e.g., T1234) . Use with $ in 
card column 72. 

or REWIND. The Write-Tape will 

be rewound before being used, 
or LIST. The Write-Tape will be 

listed by Subroutine LTAPE 
after the matrix has been 
written on the Write-Tape, 
or Anything else will be ignored. 

I The Write-Tape Number (e.g., 

21 ). 

or Blank if the matrix is not to 

be written on tape. 

Note (1) Format Type A allows any keypunch symbol. 

vormat Type I allows only integer numbers right justified 
in the field. 

As examples of tape input to Subroutine YREAD consider: 

Example 1. A matrix named AB2 with run number of RUN-46 Is to 
be read from tape number 11 into the computer and 
printed. This matrix is also to be written on tape 
number 22 that is to be initialized and identified as 
T4321. 

Example 2. A matrix named XYZ4 with run number of TKD is on tape 
number 13 twice. The first time is at location 29 
and the second time is at location 54. It is desired 
to read the second matrix. 

Figure 2 demonstrates how these two examples would be written on 
a coding form to facilitiate keypunching to cards. 




fORMA 
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Figure 2 Examples of Tape Input for Subroutine YREAD 
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Subroutine YREAOO reads a matrix from octal numbers on cards* Matrix 
[a] is then stored on NITIA in FORMA sparse notation. The matrix is printed 
side by side in both octal and decimal so that these input data are re- 
corded with the ansvers of a run. 

The primary purpose of Subroutine YREAOO is to read a matrix from punched 
cards without round-off error. The cards are punched by Subroutine YFNCHO. 
Octal representation of the matrix elements is used because it gives an 
exact replica of the binary number used by a digital computer. A decimal 
representation will not give an exact replica. The matrix on punched cards 
is to be used only as an emergency backup for the matrix written on a stor- 
age tape* 

Because of the emergency backup nature of input data to this Subroutine 
YREAOO, only cards are read. No tape re£.ding or writing options are 
available* 


The information entered on the data cards is given below. Required en- 
tries are denoted by an * symbol. Any other entry is optional. 



Card 

Columns 

Format 
Type (1) 

Entry 

First Card 

1-6 

A 

^Matrix Name. (Will appear in 


7-10 

I 

printout.) 
^Matrix Row Size. 


11-15 

I 

^Matrix Column Size. 


16-21 

A 

^Matrix Shape. 


22-69 

A 

Any remarks to further identify 

Middle Cards 

1-5 

I 

the input matrix. 

*Row Number of matrix elements 


6-10 

I 

on card. 

^Column Number of matrix elemen 


14-25 

0 

in first data field. 

^First data field with matrix 


29-40 

0 

elements. (2) 

^Second data field with matrix 


44-55 

0 

elements. (2) 

*Third data field with matrix 

Last Card 

1-10 

I 

elements. (2) 
*Ten Zeroes. 


Note (1) Format Type A allows any keypunch symbol. 

Format Type I allows only integer numbers right justified in 
the field. Format Type 0 allows only octal numbers. 

Note (2) Only non-zero elements need be entered. 
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No exanples of input are given because data would not be keypunched for 
input tc Subroutine YREADO but rather obtained from Subroutine YPNCHO. 
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Subroutine YRBVAD rearranges (revises) the rove and coluams of a aatrlx 
[a 3 * nultlplles [a 3 by a scalar alpha, and adds the result to a previously 
defined statrlx ^(0 . Ca 3 , Cz in3 , and [Z outi ara stored on NUTA and 
NUTZ In FORMA sparse notation. The revision of C^J la specified by two 
vectors. The first vector {iVECj gives the new row location of each row 
of C<^3 In CzJ • The second vector jjVEC} gives the new colunn location 
of each coluon of Ca 3 In [Z3 . The ^^3 matrix must be defined before the 
use of this subroutine. For instance, if [t] is to be originally defined 
as all zeros. Subroutine YZERO could be used. The YREVAD operation can be 
thought of In subscript notation as 

2 ^ (out) - (in) + « Oy 


idiere 


k » IVEC(i) , 
i = JVEC(J) . 

NBA is the number of rows of C^3 » and NCA is the number of columns 
of Ca3 • 

Values In |lVEC[ and {jVEcj may be positive, negative or zero. A nega- 
tive value changes the sign of the corresponding row or coluosi of [AJ in 
fZ3 . [A3zero va^ue omits the corresponding row or column of C^3 ftom 
The values are integer numbers. 

This subroutine may be called repeatedly to form Cz3 from the revision/ 
addition of several Ca 3 matrices. 

An important use of Subroutine YREVAD is to revise and add the stiffness 
matrix of a structural component to the stiffness matrix of the complete 
structure to account for the difference in coordinate systems. 


/i = 1, NRA\ 
» 1, NCAj 


EXAMPLES 


The first example to illustrate the YREVAD op ration is as follows. 
Matrix C^3 has been previously defined as 



1 . 

0 . 

0 . 

0 . 


0 . 

2 . 


0 . 

0 . 

3. 

0 . 


0 . 

0 . 

4. 

0 . 


0 . 

0 . 

5. 

6 . 


0 . 

0 . 
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Let £a 3 be defined as 



The first row of [a] Is to be added to the third row of Tzl , the second 
row of Is to be omitted from ^Z] , the third row of is to be added 

to the first row of , with the sign of each element reversed. 


Thus jlVEcl^^j - 



The first column of is to be added to the second column of with 
the sign of each element reversed, the second column of ^a 3 is to be 
added to the fifth column of • 


Thus {JVEC|2^j 

Then, assuming o -> 1.0 and placing jlVEC [and jjVEcf adjacent to [A] 
to aid In visualising tb^ revision of C^3 > ^ have 
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A second example of the use of this subroutine demonst rates 
the coordinate transformation concept which is a very important 
application of YREVAD. It is desired to transform a stiliness 
matrix from an original coordinate system (subscript 1) to a 
final coordinate system (subscript 2) as shown in tfie sketch below. 





For this example, the stiffness matrix in tlie original coordinate 
system is assumed to be 
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From inspection of the coordinate system axes in the above sketcli. 





y^ 


-^1 

J'K. 


yi 


To obtain the vector on the right hand side of the equation: 


xi of the original coordinate system is to be the second 
row of the final coordinate system with the sign reversed. 


y\ of the original coordinate system is to be the third 
row of the final coordinate system. 


Z\ of the original coordinate system is to be the first 
row of the final coordinate system with the sign reversed. 


The |lVEC| to accomplish this change is 

-2 
3 

-1 


|IVEC[, 


3x1 


Apply this jlVEC^ as the |lVEc| and jjVEc} to the original stiff- 
ness matrix to obtain the final stiffness matrix. This application 
is analogous to the change in coordinates made in the triple matrix 
product procedure. 






-2 


1. 2. A. 


6. A. -5. 

3 


2. 3. 5. 

YREVAD^ 

A. 1. -2. 

-1 


A. 5. 6. 


-5. -2. 3. 


The above YREVAD procedure may be compared witii the more conven- 
tional triple matrix product procedure below. 
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The equation for the coordinate transformation would ho 


Xl 


-1. 


^^2 

>'l 

ts 

1. 


yz 

Zl 


-1. 


Z2 


(1 


The strain energy expression with the stiffness matrix is 
given by U ® Ih. |q| or for coordinate system 1, 


U = -*5 (xj yi zj] 


1. 2. 

2 . 3 . 

4. 5. 


4. 

5. 

6. 


XI 

yi 

Z] 


( 2 ) 


Substituting Equation (1) into (2) gives 



-1. 


1. 2. 4. 


-1 . 


*2 

* h [x2 yz za] 

-1. 


2. 3. 5. 


1. 


yz 


1. 


4 . 5 . 6 . 


-1. 


-K 


thus the inner triple matrix product gives the stiffness matrix 
in coordinate system 2, i.e., 

6. 4. -5. 

4. 1. -2. 

[-5. -2. 3. 

which is the same result as that obtained from the YREVAD procedure. 

The advantages of using the YREVAD procedure over the triple 
matrix product procedure are: 


1) Less computer time; 

2) Less computer core is used; 

3) Usually easier to code the |lVLc| (lluis |jVLC|) 
than to code th.e transformation matrix. 
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DESCRIPTION OF TECHNIQUE 

Matrix Ca3 Is read from NUTA and the non-zero elenent location numbers 
are replaced according to the contents of |lVEC| and {JVEC [ . The sign 
of the non-zero elements Is changed accordingly and the revised [*S] Is 
stored on NUTl. Matrix ln3 is read from NUTZ and stored on NUT2. 
Subroutine YAABB Is then called to sum C'Xj and [Z Inj . !D)e resulting 
matrix Cz ouO -s stored on NUTZ. 



YRTAPE 


Subroutine YRTAPE reads partitions of a selected FORMA sparse 
notation matrix from a FORMA reserve tape (disc) Into the com- 
puter core and stores them sequentially on a utility tape to be 
operated on by other FORMA subroutines. The matrix to be selected 
Is Identified by the desired run number and matrix name. This 
procedure is accomplished by searching the matrix headings (See 
Subroutine YWTAPE explanation.) until a match with the desired 
run number and matrix name Is obtained, and then reading the matrix 
elements and element locations Into core and storing them on a 
utility "ape. This procedure Is repeated for each matrix partition. 
The search starts from the current position (does not rewind) of 
the tape (disk) and proceeds to the EOx (end of tape defined In 
subroutine YWTAPE explanation) . If the desired matrix was not 
found upon reaching the EOT, a rewind Is pat' jrmed and one more 
search to the EOT Is made. If the desired matrix is again not 
found, (1) an error message Is printed, (2) a listing of the 
matrix headings Is printed (See Subroutine LTAPE writeup) , rucd 
(3) transfer Is made to Subroutine ZZBOMB where the program Is 
terminated. 



YRVADl 


Subroutine IRVADl performs a special case of th3 function performed by 
FOIOfA Subroutines REVADD and YREVAD (see subroutine explanations for 
these subroutines) . Subr mtlne YRVADl rearranges (revises) the rovs 
and columns of a matrix I A] , multiplies |,A) by a scalar alpha, and adds 
the result to a previously defined matrix Z inj . Qa] is a matrix in 
FORMA dense storage. Cz In? and {Z out3 are stored on MUTZ In FORMA 
sparse notation. YRVADl may be used only tdien ^A3 , CZ Inl , and Tz outj 
are all symetrlc. Because of this symmetry, only one vector }lJVEC| is 
required for the revision of matrix fA? . 


DESCRIPTION OF TECHNIQUE 

Any non-zero elements in are stored in workspace V. Element location 
numbers are formed for the corresponding elements In workspace LV. These 
elements of Ca 3 are then stored on NUTl in FOISfA sparse notation. Matrix 
Cz InJ is read from NUTZ and stored on NUT2. Subroutine YAABB Is called 
to add sparse matrix Ca 3 to matrix fZ in3 and the resulting matrix QZ out3 
is stored on NUTZ. 

Limitation: 

The dimension size of V and . must exceed or be equal to the number of 
non-zero elements in matrix [nj . 



YRVAD2 


Subroutine v' , ad 2 reads small dense matrices and Integer vectors 
from a utility tape (disk) and uses the vectors to re''lse the 
matrix locations to form large sparse matrices. These sparse 
matrices are then summed to form one large matrix in FORMA sparse 
notation [Z]. Because the small matrices are assumed to be sym- 
metric, [Z] is also symmetric and only the lower half is formed. 
Only matrix elements having an absolute value greater than 10*^^ 
are summed in [Z]. Subroutine YRVA02 reads matrices from the 
utility tape until a blank record is found. 



YRVISl 


Subrouclne YRVISl rearranges (revises) the colusns of a matrix in 
FOniA dense storage to form matrix in FORMA sparse notation. The 
rows of fA3 are not rearranged, thus only one integer vector )JVEC| 
need be input to this subroutine. 


DESCRIPTICTi OF TECHNIQUE 

The non-zero elements of C^3 are stored in workspace V. Element location 
numbers are formed in workspace LV corresponding to the elements in V. 

Cz3 , composed of these elements and locations, is stored on NUTZ. 

Limitation: 

The dimension size of V and LV must exceed or be equal to the number of 
non-zero elements in C^3 • 
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Subroutine YSRED2 operates on the stiffness matrix [A] to form .1 
reduced stiffness matrix [R] and/or the reducing transformation 
(T). The relation between the stiffness matrix fAj, displacements 
{Xl. and applied fjrces {B} may be expressed in matrix form as 

[AJ (X) - {B) H J 

The reduction method assumes Eq [IJ to be partitioned as 

12J 


where {X^} are the displacements to be reduced out and {X^} 

are the displacements to be retained. The applied forces acting 
on the coordinates to be reduced are assumed to be zero, such that 

{B^} * {0} 

Substituting Eq [3] Into Eq [21 and expanding the upper partition^ 
will yield the reduced displacements in terms of the retained 
displacements as 

(Xj,} = - lA^2l tx^l 

Expanding Che lower partitions of Eq [2] and substituting Eq (4] 
will yield the reduced stiffness matrix as 


[31 


[41 



[RJ {X^} = {B2> 

where [K] is the reduced stiffness matrix and is expressed as 

[RJ = [A22J - [A21) [Aiil'^ [A^ 2 l 



15J 


[ 6 J 


[T] {X^} 


[7] 
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where 


m - 


Also 

(R) - lAj m (91 

DE SCRIPTION OF T ECHNIQUE 

This subroutine uses Gauss reduction partially completed to form 
matrix [R) and [T] from stiffness matrix [A]. As an example of 
the method, consider three simultaneous equations of the following 
form 


•lAjil-' lAjjI 

m 


[a] 


®11 *1 ®12 *2 ®13 ’‘3 “ 

“21 *1 ^ “22 *2 ^ ®23 *3 “ *»2 

*1 ^ ^32 ^2 ^33 ^^3 = ^ 


tlOJ 


These equations may be written in matrix form as 


®11 ®12 ®13 


®21 ®22 ®23 


®31 ®32 ®33 



rbp 




■2 


'3J 


[111 


Solve the first Equation for x^ as 

®12 ®13 . 

1 2 3 a^. 


(121 


Substituting £q [12] into the second and third equations in Eq [11] 
and divide the first by results in 

- 1 

* A 

1 a^2 ®13 

it ★ 

“ “22 =23 

it it 

0 a^2 833 



[131 



where 


* 

“12 

®12 " 

“11 

* 

“13 

^13 * 

“11 


®22 “ ®22 


^12 ^21 
®11 


®23 “ ®23 


^13 ^21 


“32 “ “32 


“l2 “31 
®11 


“33 “ “33 

K* ‘’1 

b - — 

^ “11 


“l3 “31 

“11 


b2 = b2 - 


\ “21 
“11 


b3 = b3 - 


^ “3 1 
“11 
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114! 

U5I 

1161 

117] 

1181 

1191 

120 ) 

1211 

122 ) 


Solve the second equation for 

“*3 . _^2 

^2 ~ * *3 ‘ 

“22 “22 


which will yield 


1231 


Substitute £q [23) into the third equation in Eq [13) and divide 
the second equation by ^ 22 ’ results in 


H- 1 

* 

“12 

A 

“13 


r > 


r 

* ^ 

0 

1 

** 

“23 

< 

^2 


1 

‘>2 

0 

0 

** 

“33 



1 

irk 

bJ 


124) 
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where 


** 

*23 


■23 

V 

*22 


A* 

‘33 


AA 


A A 


*23 32 

®33 * 

®22 


22 


AA * 

63 = 63 


A A 
**2 ®32 


*22 


125J 


126J 


(27J 


(281 


The reduced stilfness matrix has been formed and is contained as the 

AA 



1291 

(301 




[A12I 


(A21I 


ii 

*12 

‘21 

*22 


■13 


23 


(a 3 i 332! 



(31] 

(32) 

(331 

(34) 

(351 

(361 
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Substituting Eq 131] through [3A] into Eq [6] results in a reduced 
stiffness matrix of the form 


. 21 ®23 31 12 13 *32 13 22 31 

IRJ ■ a__ + 

33 a. . a,, - a.^ a. 


“ll *23 *32 


11 22 


‘12 “21 


137] 


Equation [37] is Identical to the result obtained by expanding 
Eq [26]. Thus, Gauss reduction partially completed yields the 
reduced stiffness matrix. 

The reducing transformation may also be obtained using Gauss re- 
duction If additional operations are performed. From Eq [24], 
solve the second equation for 


A* icif 

*2 " ' “23 ‘3 * ‘'2 

Substitute Eq [38] into the first equation in Eq [24], which will 
yield 




r ^ 


r **;ti 

1 0 a^3 




‘’i 

** 

i 


r< 


0 1 a23 

< 

*2 


Alt 




AA 

0 0 




b. 

33 


3 

^ J 


L ^ J 


where 


Hifit it * ** 

*13 " *13 " *12 *23 

* 5*C ** 

\ ’ h ■ *12 '’2 


[38] 


[39] 


[40] 

[41] 


Inspection of Eq (39] .shows that we have formed a unity matrix 
partition wliere the row-columns were reduced* The afi!* ^ind the 

a 2 j elements contain the necessary Information to form the reduc- 
ing tccinsformation. To show this, substitute Eq [31] and [32] into 
Eq [8J to yield 



[TJ = 


®21 

.^2JL 

- “13 

“22 


^22 

- 

“12 

"*12 

:‘l3 



“ll 

"'*22 

<N 

1 

“12 
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1^21 


The second and third rows of Eq [42] are equal to the negative uf 

ititit 

elements and a^^ in Eq [39]. Thus, Gauss reduction also 

yields the reducing transformation. 


DESCRIPTION OF TECHNIQUE 

The stiffness matrix ^a 3 is read from NUTA in FORHA sparse notation, 
converted to banded notation, and stored on NUTl. This banded nota- 
tion stores groups of colunns of the upper half of synnetric matrix 
[a 3 . The reducing procedure is similar to the decomposition proce- 
dure of subroutines YDC0M3 and YDCM3A. ^le reduced stiffness matrix 
fRj is stored on NUTR in FORMA sparse notation. If input argument 
"in" is equal to one, the reducing transformation is formed and 
stored on MUTT in rORMA sparse notation. 



YSTOD 


Subroutine YDTOS reads a FORMA sparse matrix from the peripheral 
device on which the matrix Is written and stores It In a dimensional 
work space In core. The work space is zeroed before the nonzero 
terms of the sparse matrix are Inserted into the work space. 

DESCRIPTION OF TECHNIQUE 

Dense matrix [E] Is formed from the nonzero terms of FORMA sparse 
matrix [A]. After all of [Z] Is zeroed, each partition of [A] Is 
sequentially read Into work vectors {V} and (LV). For each non- 
zero element location, LV(k), 1 and J are found: 

1 « LV(k)/100000 
J “ LV(k) -100000 * 1 

then 

z^j = V(k) for all elements of [A]. 



YSYMUH 


Subrouclne YSYMUH symmetrizes a square FORMA sparse matrix [A]. 
Existing elements above the diagonal are set to zero and new ele- 
ments are added above the diagonal to reflect elements below the 
diagonal. The matrix Is then searched and any zero elements are 
removed. Finally, the elements and their locations are ordered. 
This result Is [zj where: 


Z * A 

ij iJ 

i 1 J 

and 



i < J . 



YTRANS 


Subroutine YTRANS calculates the transpose (Interchange of rows 
and columns) of a matrix In FORMA sparse notation. If 

Is the matrix to be transposed, then the result Is 

^^^NCAxNRA 

where 

/I - 1, NRA\ 

*jl “ ®lj \J = 1. nca/ 

NRA is the number of rows of [A], and NCA Is the number of columns 
of [A]. 

DESCRIPTION OF TECHNIQUE 

The location of each a (100000 * 1 + j) Is transposed to 
100000 * J + 1 and the matrix is reordered. 



YUNITY 


Subroutine YUNITY generates a square matrix [Z*] with diagonal elements 
equal to one and all o£f*dlagonal elements equal to zero. That Is, 


-1. (1 “ J) 

z^j - 0. it * J) 

In matrix notation, 


[Z] 


NxN 


1 . 


0 . 


0 . 


1 . 


where N la the sire of [Z] (square) . A synonym for the unity matrix 
la the Identity matrix, thus the usual designation as • 

A matrix Is unaltered when multiplied by the unity matrix and the pro- 
cess is connutatlve. In matrix notation. 


H W • W W - M • 

nE8r.Mmon op TEcamouE 

N ones are stored in workspace V an-l N diagonal element location numbers 
are stored In workspace LV. , composed of these ones and location 

numbers Is stored on NUTA. 



YWRITE 


Subroutine WRITE vrltes a FORMA bparse notation matrix of real 
numbers (a FORTRAN term for numbers with a decimal point) on paper. 
A group of up to ten consecutive elements from a row of the matrix 
Is printed on each line. If all of the elements of a group are 
zero, printing of this line is suppressed. 

Each matrix printed begias on a new page. On each page of print- 
out is the page heading given by Subroutine PAGEHD, the name 
of the matrix, the row sli:e and column size of the matrix, the 
number on nonzero elements In the matrix, and the number of parti- 
tions In the matrix. 

This Is followed by the matrix data. On any line of matrix data 
the first integer number is the row number of the matrix elements 
on that line. The second Integer number Is the coluiun number of 
the matrix element In the first data field. The next group of 
real numbers (up to ten) ^re the values of the matrix elements. 
This group of matrix elements Is given In consecutive column 
order. 
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Subroutine YWTAPE writes FORMA sparse matrix data at the end of 
existing written matrix data on a FORMA tape (disk is preferred^ 
see below). Each sec of matrix data consists of two logical 
records. The first record contains the matrix heading (tape 
identif ication» location number^ run number of rows of matrix^ 
number of columns of matrix* date* the acronym "spare*** the 
number of nonzeroes in ^ partition* the number of the partition* 
and the number of partitions in the matrix* The second record 
consists of the matrix elements of the partition and the element 
locations* Subroutine YWTAPE is compatible with FORMA subroutine 
WTAPE and dense and sparse matrices may be stored on the same 
FORMA reserve tape* 

A schematic representation of the tapr (disk) is given by the fol- 
lowing sketch* 


Beginning 

of 

tape (disk) 


H* 

E* 

H., 

E, 


H 

E 

EOT 

s 

1 

1 

2 

i. 


n 

n 


2 


where 


= Matrix h .ading of the i^*^ written matrix partition* 

= Matrix elements of the i^^ written matrix partition* 

EOT = End of Tape. Data written by Subroutine WTAPE, YWTAPE* or 
INTAPE that all FORMA tape subroutines recognize as being 
the end of written data. 

Each vertical line is au end of logical rec'^rd puw on by computer 
system *s routines. The tape Is written xn binary form as ^'pposed 
to binary coded decimal (BCD) form. 

To find the end of written m;»trix data* a search is made of the 
matrix headings until the EOT is found. For this reason* a 
"new" cape (disk) muse be initialized with Subroutine lOTAPE so 
that the tape (d^sk) contains an EOT. A **new" tape (dxsk) is 
defluzJ to be a tape (disk) for which it is desired to start 
writing matrix data at the front of the tape (disk). Thus, a 
"new" tape (disk) could be om vith obsolete FORMA matrix data 
on it as well as one that ha- ivar been written on by whe 
FORMA system. When the EOT found* a backspace operation is 
done over tha EOT* and then the current marrix heading, current 
matrix elements, and a nev EOT is written. 
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A disk is preferred to a tape because the physical separation of 
the read and write heads on most tape drives may cause tape 
tolerance problems; thus back-spacing over the EOT is usually not 
successful. Instead of ending up positioned in front of the EOT, 
the write head is often positioned in front of the previous matrix 
elements (E^ is the ai ^ve sketch). The current matrix heading 

will be written over previous matrix elements. This causes 
problems later when trying to read th'^ records written on the tape. 
To alleviate this problem, it is strongly reconinended that all 
FORMA tape subroutines use an intermediate device such as a disk. 

At the start of a computer run, the existing tape should be copied 
on the disk by using computer control cards. Likewise, at the 
end of the run, the disk should be copied back on tape by using 
computer control cards. 



YZERLH 


Subroutine YZERLH elloinates the non-zero elements below 
of a matrix C^3 stored on NUIA in FURM& sparse notation. 


a^j - 0. (i > J) 


EXAMPLE 


If [Aj is input 


to Subroutine YZERLH as 



2 . 3 . 
5. 6. 

8. 9. 


the matrix output from this subroutine will be 



3. 

6 . 

9- 


the diagonal 
That is. 


Matrix Qa 3 is read into workspaces V and LV. The non-zero elemertt 
locations are searched for element locations below the diagonal. These 
locations and the corresponding elements are eliminated. 



YZERim 


Subroutine YZERUH elioinatee the non-zero elenents above the diagonal of 
a matrix Ca] stored on NUIA In Forma sparse notation. That Is, 

=0. (1 < j) 

EXAMPLE 

If CaJ la Input CO Subroutine YZERUH as 



the matrix output from this subroutine will be 



Matrix Is read Into trorkspaces V and LV. The non-z..ro element 
locations are searched for element locations above the diagonal. These 
locations and the corresponding elements are eliminated. 



Y2ER0 


Sub rout Luo YZCRO genera tes a matrix with ea*h equal to 

zero. That is, 


2 . . - 0 . 

J-J 


c 


I. NK\ 

1 , Nc; 


In matrix notation, 

b. 0 . ... 0. 

0 . 0 . 

0 . . . . 0 . 

where NR is the number of rows of (Zj, and NC is the auiaher of 
columns of [Z]. 


[y^] 


NRxNC 


This subroutine is useful in setting a matrix array to zero 
before perfonaing subsequent operations such as matrix assembly 
(jlASSEM) or rcvision/addi tion of one matrix Into anothe- (YREVAD). 

Matrix stored on NUTA in FORMA sparse notation. Because 

matrix elements equal to zero are not stored in thl. notation, only 
matrix headers are stored on NUTA. 



